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Abstract 

The Planck High Frequency Instrument (HFI) is designed to measure the temperature and polarization anisotropies of the Cosmic 
Microwave Background and galactic foregrounds in six wide bands centered at 100, 143, 217, 353, 545 and 857 GHz at an angular 
resolution of 10' (100 GHz), 7' (143 GHz), and 5' (217 GHz and higher). HFI has been operating flawlessly since launch on 14 May 
2009. The bolometers cooled to 100 mK as planned. The settings of the readout electronics, such as the bolometer bias current, 
that optimize HFI's noise performance on orbit are nearly the same as the ones chosen during ground testing. Observations of 
Mars, Jupiter, and Saturn verified both the optical system and the time response of the detection chains. The optical beams are 
close to predictions from physical optics modeling. The time response of the detection chains is close to pre-launch measurements. 
I The detectors suffer from an unexpected high flux of cosmic rays related to low solar activity. Due to the redundancy of Planck's 
observation strategy, the removal of a few percent of data contaminated by glitches does not significantly affect the sensitivity. The 
I cosmic rays heat up the bolometer plate and the modulation on periods of days to months of the heat load creates a common drift of 
I all bolometer signals which do not affect the scientific capabilities. Only the high energy cosmic ray showers induce inhomogeneous 
• • heating which is a probable source of low frequency noise. The removal of systematic effects in the time ordered data provides 
^ a signal with an average level of noise less than 70% of our goal values in the 0.6-2.5 Hz range. This is slightly higher than the 
pre-launch measurements but better than predicted in the early phases of the project. This is attributed to the low level of photon 
noise resulting from an optimized optical and thermal design. 
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1. Introduction 

Planck 1 (Tauber et al. 2010a; Planck Collaboration 2011a) 

* Corresponding author: J.-M. Lamarre, jean-michel. 
lamarre@obspm.fr. 

1 Planck (http://www.esa.int/PZancA:) is a project of the 
European Space Agency (ESA) with instruments provided by 
two scientific consortia funded by ESA member states (in par- 



is the third-generation space mission to measure the 
anisotropy of the cosmic microwave background (CMB). 
It observes the sky in nine frequency bands covering 30- 



ticular the lead countries France and Italy), with contributions 
from NASA (USA) and telescope reflectors provided by a collab- 
oration between ESA and a scientific consortium led and funded 
by Denmark. 
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857 GHz with high sensitivity and angular resolution from 
31' to 5'. The Low Frequency Instrument (LFI; Mandolesi 
et al. 2010; Bersanelli et al. 2010; Mennella et al. 2011) cov- 
ers the 30, 44, and 70 GHz bands with amplifiers cooled to 
20 K. The High Frequency Instrument (HFI; Lamarre et al. 
2010; Planck HFI Core Team 2011a) covers the 100, 143, 
217, 353, 545, and 857 GHz bands with bolometers cooled 
to 0.1 K. Polarization is measured in all but the highest two 
bands (Leahy et al. 2010; Rosset et al. 2010). A combina- 
tion of radiative cooling and three mechanical coolers pro- 
vides the temperatures needed for the detectors and optics 
(Planck Collaboration 2011b). Two data processing centres 
(DPCs) check and calibrate the data and make maps of the 
sky (Planck HFI Core Team 2011b; Zacchei et al. 2011). 
Planck's sensitivity, angular resolution, and frequency cov- 
erage make it a powerful instrument for galactic and ex- 
tragalactic astrophysics as well as cosmology. Early astro- 
physics results are given in Planck Collaboration (2011h-x). 

The goal of this paper is to describe the in-flight perfor- 
mance of the HFI in space and after the challenging launch 
conditions. It does not attempt to duplicate the content of 
the Planck pre-launch status papers (Lamarre et al. 2010; 
Pajot et al. 2010), but rather presents the operational sta- 
tus from an instrumental viewpoint. These results propa- 
gate to scientific products through the data processing re- 
ported in the companion paper (Planck HFI Core Team 
2011b) which describes the instrumental properties as they 
appear in the maps used by the "Planck early results" com- 
panion papers. This paper focuses on the ability of the HFI 
to measure intensity without any description of its perfor- 
mance in measuring polarization, which will be reported 
later. 

Section 2 summarizes the instrument design. Section 3 
focuses on early in-flight operations, the verification phase 
and the setting of the parameters that have to be tuned in 
flight. Section 4 addresses the measurement of the beams on 
planets and the disentangling of time response effects from 
the beam shape. It also presents the best current knowledge 
of the physical beams resulting from this work. The effec- 
tive beam obtained after data processing are to be found 
in Planck HFI Core Team (2011b). Sections 5, 6 and 7 are 
dedicated to noise, systematic effects and instrument sta- 
bility respectively. A summary of the HFI in-flight perfor- 
mance and a comparison with pre-launch expectations are 
presented in section 8. 

2. The HFI instrument 

2.1. Design 

The High Frequency Instrument (HFI) was proposed to 
ESA in response to the announcement of opportunity for 
instruments for the Planck mission in 1995. It is designed to 
measure the sky in six bands (Tab. 1) with bolometer sen- 
sitivity close to the fundamental limit set by photon noise. 
The lower four frequency bands include the measurement 
of the polarization. This sensitivity is obtained through a 
combination of technological breakthroughs in each of the 
critical components needed for bolometric detection: 

— Spider web bolometers (Bock et al. 1995; Holmes et al. 
2008) and polarization sensitive bolometers (Jones et al. 
2003) which can reach the photon noise limit with suf- 
ficient bandwidth to enable scanning great circles on 



the sky at roughly 1 rpm. They offer a very low cross- 
section to cosmic rays that proves to be essential in this 
environment and with this sensitivity. 

— A space qualified 100 mK dilution cooler (Benoit et al. 
1997) associated with a high precision temperature con- 
trol system. 

— An active cooler for 4K (Bradshaw & Orlowska 1997) 
using vibration controlled mechanical compressors to 
prevent excessive warming of the 100 mK stage and min- 
imize parasitic effects on bolometers. 

— AC biased readout electronics that extend high sensi- 
tivity to very slow signals (Gaertner et al. 1997). 

— A thermo-optical design consisting, for each optical 
channel, of three corrugated horns and a set of compact 
reflective filters and lenses at cryogenic temperatures 
(Church et al. 1996). These include high throughput 
(multimoded) corrugated horns for the 545 and 857 GHz 
channels (Murphy et al. 2002). 

The angular resolution was chosen to extend the mea- 
surement of the small scale features in the CMB, while 
keeping the level of stray light to extremely low levels. At 
the same time, at this sensitivity, the measurement and 
removal of foregrounds requires a large number of bands 
extending on both sides of the foreground minimum. This 
is achieved with the six bands of the HFI (Table 1) and 
the three bands of the Low Frequency Instrument (LFI; 
Mennella et al. 2011). 

The instrument uses a ~ 20 K sorption cooler com- 
mon to the HFI and the LFI (Planck Collaboration 2011b; 
Bhandari et al. 2000, 2004). The HFI focal plane unit 
(FPU) is integrated inside the mechanical structure of the 
LFI, on axis of the focal plane of a common telescope 
(Tauber et al. 2010a). 

The ability to achieve background limited sensitivity 
was demonstrated by the ARCHEOPS balloon-borne ex- 
periment (Benoit et al. 2003a,b) , an adaptation of the HFI 
designed for operation in the environment of a stratospheric 
balloon. Similarly, the method of polarimetry employed by 
the HFI was demonstrated by the Boomerang experiment 
(Montroy et al. 2006; Piacentini et al. 2006; Jones et al. 
2006) . The HFI itself was extensively tested on the ground 
during the calibration campaigns (Pajot et al. 2010) at IAS 
in Orsay and CSL at Liege. However, the fully integrated 
instrument was never characterized in an operational envi- 
ronment like that of the second Earth-Sun Lagrange point 
(L2) . In addition to thermal and gravitational environmen- 
tal conditions, the spectrum and flux of cosmic rays at L2 
is vastly different than that during the pre-flight testing. 
Finally, due to the operational constraints of the cryogenic 
receiver, the end to end optical assembly could not be tested 
on the ground with the focal plane instruments. 

The instrument design and development are described 
in Lamarre et al. (2010). The calibration of the instrument 
is described in Pajot et al. (2010). The overall thermal and 
cryogenic design and the Planck payload performance are 
critical aspects of the mission. Detailed system-level aspects 
are described in Planck Collaboration (2011a) and Planck 
Collaboration (2011b). 

2.2. Spectral transmission 

The spectral calibration is described in Pajot et al. (2010) 
and consists of pre-launch data, in the passband and 
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Figure 1. HFI spectral transmission 



around, combined with component level data to determine 
the out of band rejection over an extended frequency range 
(radio-UV). Analysis of the in-flight data shows that the 
contribution of CO rotational transitions to the HFI mea- 
surements is important. An evaluation of this contribution 
for the J = 1 -> (100 and 143 GHz), J = 2 -> 1 (217 
GHz) and J = 3 — > 2 (353 GHz) transitions of CO is pre- 
sented in Planck HFI Core Team (2011b). 



3. Early HFI operation 

3.1. HFI Cool down and cryogenic operating point 

The Planck satellite cooldown is described in Planck 
Collaboration (2011b). 

The first two weeks after launch were used for passive 
outgassing, which ended on 2 June 2009. During this pe- 
riod, gas was circulated through the 4 He-JT cooler and the 
dilution cooler to prevent clogging by condensable gases. 
The sorption cooler thermal interface with HFI reached a 
temperature of 17.2 K on 13 June. The 4 He-JT cooler was 
only operated at its nominal stroke amplitude of 3.5 mm 
on 24 June to leave time for the LFI to carry out a spe- 
cific calibration with their reference loads around 20 K. The 
operating temperature was reached on 27 June, with the 
thermal interface with the focal plane unit at 4.37 K. 

The dilution cooler cold head reached 93 mK on 3 July 
2009. Taking into account the specific LFI calibration re- 
quirement that slowed down the cooldown, the system be- 
haved as expected within a few days, according to the ther- 
mal models adjusted to the full system cryogenic tests in 
the summer 2008 at CSL (Liege). 

The regulated operating temperature point of the 4K 
stage was set at 4.8 K for the 4K feed horns on the FPU. 
The other stages were set to 1.395 K for the so called 
1.4 K stage, 100.4 mK for the regulated dilution plate, and 
103 mK for the bolometer regulated plate. 



These numbers were very close to the planned operat- 
ing point. As the whole system worked nominally, margins 
on the cooling chain for interface temperatures and heat 
lift are large. The Planck active cooling chain was one of 
the great technological challenges of this mission and is 
fully successful. A full description of the performance of 
the cryogenic chain and its system aspects can be found 
in Planck Collaboration (2011b). The parameters of the 
operating points of the 4K, 1.4 K and 100 mK stages are 
summarised in Table 2. 

The temperature stability of the regulated stages has 
a direct impact on the scientific performance of the 
HFI. These stabilities are discussed in detail in Planck 
Collaboration (2011b). Their impact on the power received 
by the detectors is given in Sect. 3.3.1. 

3.2. Calibration and performance verification phase 
3.2.1. Overview 

The calibration and performance verification (CPV) phase 
of the HFI consisted of activities during the initial cooldown 
to 100 mK and during a period of about six weeks before 
the start of the survey. The cooldown phase is summarized 
in Sect. 3.1. The pre-launch value of the 4 He-JT cooler op- 
erating frequency was used (see Sect. 3.2.2). Activities re- 
lated to the optimization of the detection chain settings 
were performed first during the cooldown of the JFET am- 
plifiers, and again when the bolometers were at their oper- 
ating temperature. Most of the operating conditions were 
pre-determined during the ground calibration. The main 
unknown was the in-flight background on the detectors. 
The detection chain settings are presented in Sect. 3.2.3. 
Other CPV activities performed are: 

— determination of the detection chain time response un- 
der the flight background 

— determination of the detection chain channel-to-channcl 
crosstalk under the flight background 

— characterization of the bolometer response to the 4K 
and 1.4 K optical stages, and to the bolometer plate 
temperature variations 

— checking the immunity of the instrument to the satellite 
transponder 

— optimization of the numerical compression parameters 
for the actual sky signal and high energy particle glitch 
rate 

— various ring-to-ring slew angles (1.'7, 2.'0 [nominal], 2.'5) 

— checking the effect of the scan angle with respect to the 
Sun 

— checking the effect of the satellite spin rate around its 
nominal value of 1 rpm. 

On 5 August 2009, an unexpected shutdown of the 4 He- 
JT cooler was triggered by its current regulator unit (CRU). 
Despite investigations into this event, its origin is still un- 
explained. A procedure for a quick restart was developed 
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Table 2. Main operation and interface parameters of the cooling chain 

Interface Sorption cooler- 4 He-JT cooler (4K gas pre-cooling temperature) 17.2 K 

Interface 4 He-JT cooler-dilution cooler (dilution gas pre-cooling temperature) 4.37 K 

Interface 1.4 K cooler-dilution gas precooling 1.34 K 

Temperature of dilution plate (after regulation) 100.4 mK 

Temperature of bolometer plate (after regulation) 103 mK 

Temperature of 1.4 K plate (after regulation) 1.395 K 

Temperature of 4K plate (after regulation) 4.80 K 

Dilution plate PID power 24.3-30.7 nW 

Bolometer plate PID power 5.1-7.4 nW 

1.4 K PID power ~ 270 /iW 

4K PID power 1.7 mW 

4 He-JT cooler stroke amplitude 3450 (an 

Dilution cooler 4 He flow rate 16.19-16.65 /imole/s 

Dilution cooler 3 He flow rate 5.92-6.00 /imole/s 

Present survey life time (started 6 August 2009) 29.4 months 



and implemented in case the problem recurred, but it has 
not. Six days were required to re-cool the instrument to its 
operating point. The two- week first light survey (FLS) fol- 
lowed this recovery, starting on 15 August 2009. The FLS 
allowed assessment of the quality of the instrument settings, 
readiness of the data processing chain, and satellite scan- 
ning before the start of science operations. The complete 
instrument and satellite settings were validated and kept, 
and science operations began. All activities performed dur- 
ing the CPV phase confirmed the prc-launch estimates of 
the instrument settings and operating mode. We will detail 
in the following paragraphs the most significant ones. 

3.2.2. 4 He-JT cooler operating frequency setting 

The 4 He-JT cooler operating frequency was set to the 
nominal value of 40.08 Hz determined during ground tests. 
Once the cryochain stabilized, the in-flight behaviour of the 
cooler was found to be very similar to that observed during 
ground tests. The lines observed in the signal due to known 
electromagnetic interference (EMI) from the 4 He-JT cooler 
drive electronics have the same very narrow width. The 
long term evolution of the 4 He-JT cooler parasitic lines is 
discussed in Sect. 6. 

3.2.3. Detection chain parameters setting 

The JFET preamplifiers are operated at the temperature 
which minimizes their noise. This setting was checked when 
the bolometers were still warm (above 100 K) during the 
cooldown, since the bolometer Johnson noise was then 
much lower than the JFET noise. Optimum noise perfor- 
mance of the JFETs was found close at 130 K, in agreement 
with the ground calibration. 

After ground calibration, the only parameters of the 
REU remaining to optimize in-flight were the bolometer 
bias current and the phase of the lock-in detection, which 
slightly depends on the bolometer impedance. Fig. 2 shows 
the bolometer responses for a set of bias current values 
measured while Planck was scanning the sky. For this se- 
quence, the satellite rotation axis was fixed. For each bias 
value, the total detection chain noise was computed after 
subtraction of the sky signal. Ground measurements have 
shown that the minimum NEP and the maximum respon- 
sivity bias currents differ by less than 1%. Because of its 



higher signal-to-noise ratio, we use the responsivity to opti- 
mize the bias currents (Catalano et al. 2010). The optimum 
in-flight bias current values correspond to the pre-launch 
estimates within 5%. Therefore the pre-launch settings, for 
which extensive ground characterizations were performed, 
were kept (Fig. 2). In a similar way, the lock-in phase was 
explored and optimized, and again the pre-launch settings 
were kept. 

The optical background power on the bolometers is on 
the low end of our rather conservative range of predictions, 
even lower than expected from the ground measurements. 
This is attributed to a low telescope temperature and no 
detectable contamination of the telescope surface by dust 
during launch. This should result in a level of photon noise 
lower than initially expected and an improved sensitivity. 

3.2.4. Numerical data-compression tuning 

The output of the readout electronic unit (REU) consists of 
one number for each of the 72 science channels for each half- 
period of modulation (Lamarre et al. 2010). This number, 
Sreui i s t ne exact sum of the 40 16-bit ADC signal values 
obtained within the given half-period. The data processor 
unit (DPU) performs a lossy quantization of Sreu- 

First, 254 Sreu values corresponding to about 1.4 s of 
observation for each detector, covering a strip of sky about 
8° long, are processed. These 254 values are called a com- 
pression slice. The mean < Sreu > of the data within each 
compression slice is computed, and data are demodulated 
using this mean: 

Sdemod.i = (Sreu,!— < Sreu >) * ( — 1)* (1) 

where 1 < i < 254 is the running index within the com- 
pression slice. 

Then the mean < Sdcmod > of the demodulated data 
•Sdemod.i i s computed and subtracted. The resulting data 
slice is quantized according to a step Q fixed per detector: 

Sbpu,« = round ((S demod ,i- < Sdcmod >)/Q) (2) 

This is the lossy part of the algorithm: the required com- 
pression factor, obtained through the tuning of the quan- 
tization step Q, adds some extra noise to the data. For 
a/Q = 2, where a is the standard deviation of Gaussian 
white noise, the quantization adds 1% to the noise (Pajot 
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Bios [Digits] Bios [Digits] Bios [Digits] 

Figure 2. Optimization of the bolometer bias currents. Vertical lines indicate the final bias value setting. These values 
are shifted with respect to the maximum because a dynamic response correction has been taken into account. A bias 
value of 100 digits corresponds approximately to O.f nA. 



et al. 2010; Pratt 1978). In flight, the value of a was deter- 
mined at the end of the CPV phase after subtraction of the 
signal from the timeline. 

The two means < Sreu > and < Sderood > computed as 
32-bit words are sent through the telemetry, together with 
the Sdpu,! values. A variable length encoding of the Sdpu,j 
values is performed on board, and the inverse decoding is 
applied on ground. This provides a lossless transmission of 
the quantized values. A load limitation mechanism inhibits 
the data transmission, first at the compression slice level 
(compression errors) , and second at the ring level (Lamarre 
et al. 2010). 

For a given Q value, the load on each channel depends 
on the dynamic range of the signal above the level of the 
noise. This dynamic range is largest for the high frequency 
bolometers because of the galactic signal. The large rate 
of glitches due to high energy particle interactions also 
contributes to the load of each channel. Optimal use of 
the bandpass available for the downlink (75kbs _1 average 
for HFI science) was obtained by using initially a value 
of Q — cr/2.5 for all bolometer signals, therefore includ- 
ing a margin with respect to the requirement of u/Q = 2. 
The load on each HFI channel is shown and compared to 
simulated data in Fig. 3. The increase of signal gradients 



while scanning through the galactic center in September 
2009 triggered the load limitation mechanism (compres- 
sion error) and up to 80000 samples were lost for each of 
the 857 GHz band bolometers. Therefore a new value of 
Q — a/2 was set for those bolometers from 21 December 
2009 onward, reducing the number of samples lost to less 
than 200 during the following scan through the galactic cen- 
ter in March 2010. An illustration of a compression error 
loss is shown in Fig. 4. Thanks to the redundancy of the 
Planck scan strategy and the irregular distribution of the 
few remaining compression errors, no pixels are missing in 
the maps of the high signal-to-noise ratio galactic center 
regions. Periodic checks of the noise value a are done for 
each channel, but no deviation requiring a change in the 
quantization step Q has been encountered so far. 



3.2.5. Instrument readiness at the end of the CPV phase 

The overall readiness of the instrument was assessed during 
the FLS. This end-to-end test was completely successful, 
from both the instrument setting and the satellite scanning 
points of view. The part of the sky covered during the FLS 
was included in the first all sky survey. 
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Figure 3. Load measured for each HFI channel on 16 July 
2010. Simulated data for the same patch of the sky are 
shown for bolometers. Channels #54 and higher corre- 
spond to the fine thermometers on the optical stages of 
the instruments, plus a fixed resistor (#60) and a capaci- 
tor (#61) on the bolometer plate. 
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Figure 4. Example of loss in one compression slice of data 
on bolometer 857-1. Note the large signal-to-noise ratio 
while scanning through the galactic center. 



3.3. Response 

3.3.1. Variation of the signal with background and with the 
bolometer plate temperature 

The optical background on the bolometers originates from 
the sky, the telescope, and from the HFI itself. The oper- 
ating point of the bolometers is constrained by this total 
optical background, and the fluctuations of this background 
have a direct impact on the stability of the HFI measure- 
ments. 

The power spectral density of each contribution to the 
background is compared to 30% of the total noise measured 
in-flight (NEPi column of Table 6). This specification cor- 
responds to a quadratic contribution smaller than 5% on 
the total noise. 

The in-flight temperature stability of the HFI cryogenic 



Table 3. Relative response deviation (in %) from linearity 
for the CMB dipole, the galactic center (GC) and planets. 
Saturation (Sat.) occurs for the Jupiter measurements at 
high frequency. 
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stages is discussed in Planck Collaboration (2011b). The 
optical coupling of the HFI bolometers to each cryogenic 
stage is shown in the left panels of Figs. 5 and 6 and in 
Fig. 7. (The fact that the 100 mK couplings all agree with 
pre-launch measurements shows that no bolometers were 
damaged during launch.) These couplings are used to cal- 
culate the effect of the fluctuations of each cryogenic stage 
on the bolometer signals. The right panels of Figs. 5 and 
6 show the power spectral density (PSD) of the respective 
thermometers scaled by the optical coupling factors for the 
most extreme bolometers. The scaled PSDs of the thermal 
fluctuations of the 4K and 1.4 K stages are below the line 
corresponding to 30% of the total noise of the correspond- 
ing bolometer for all frequencies above the spacecraft spin 
frequency. 

The bolometer plate thermometers have a large cos- 
mic particle hit rate (Planck Collaboration 2011b) be- 
cause of the large size of their sensors compared to that 
of the bolometers. Cosmic ray hits detection and removal 
do not allow us to reach the thermometer nominal sensitiv- 
ity, therefore they cannot be used to remove the effect of 
bolometer plate temperature fluctuations on the bolometer 
signal. Instead, the data processing pipeline (Planck HFI 
Core Team 2011b) uses blind bolometers located on the 
bolometer plate. The bolometer noise components are dis- 
cussed in Sect. 6. 



3.3.2. Linearity 

The way a bolometer transforms absorbed optical power 
into a voltage is not a linear process because both the con- 
ductance between the bolometer and the heat sink, and 
the bolometer impedance have a non-linear dependence on 
the temperature (see e.g. Catalano (2008); Sudiwala et al. 
(2000)). 

The characterization of the linearity of the HFI detec- 
tors has a direct impact on the calibration of the instru- 
ment: strong non-linearity takes place during the galaxy 
crossing for high frequency bolometers and during planet 
crossings. An accurate absolute calibration is also neces- 
sary for the CMB dipole. Finally, the energy scale of large 
glitches can be corrected. The static response has been char- 
acterized during ground calibration showing a small devia- 
tion from linearity around a tenth of one percent for fainter 
sources (few hundreds of attowatts) and around a few per- 
cent for brighter sources like planets (Pajot et al. 2010). 
The static response measured during the CPV phase agrees 
with the ground estimate to better than 1%. Nevertheless, 
the use of the static non-linearity determination does not 
represent the true bolometric non-linear behaviour when 
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Figure 5. Left: coupling coefficients of the 4K stage. Right: scaled power spectral density (PSD) of the 4K stage thermal 
fluctuations for the 100- la and 353-5a-7a bolometers. 
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Figure 6. Left: coupling coefficients of the 1.4K stage. The thermal emission in high frequency bands becomes too small 
to be measured. Right: scaled PSD of the 1.4 K stage thermal fluctuations for the 100-la and 353-5a-7a bolometers. 
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times the gain (depending on the amplitude of the signal) 
and convolving it with the temporal transfer function nor- 
malized to 1 at the lowest frequency, is valid in the case of 
small signals. However this is not the case for bright point 
sources for which the estimate of non-linearity using the 
static response may be incorrect by up to 40% in the ex- 
treme case of Jupiter. For these sources we use a model 
to correct the static results. The use of fainter planets like 
Mars to characterize the beams minimizes this effect. 

Table 3 gives the deviation from linearity for various 
sources at the center of the beam for the bolometers at 
each frequency. 



3.4. Electrical crosstalk on HFI detectors 



Figure 7. Bolometer signal coupling coefficients to the 
lOOmK bolometer plate. 



scanning through bright point sources like planets. The lin- 
earization of the response done by multiplying the signal 



The electrical coupling of the signal of one bolometer into 
the readout chain of another, or electrical crosstalk, was 
measured to be less than — 60 dB for all pairs of chan- 
nels during ground-based tests (Pajot et al. 2010). We per- 
formed two tests in flight to verify this result, described 
below. 
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Figure 8. Top: Electrical crosstalk matrix Cy 54x54 for all 
bolometers, coefficients are in dB. Bottom: Distribution of 
electrical crosstalk coefficients in dB. 



3.4.1. CPV crosstalk measurements 

During the CPV phase, we switched off each readout chan- 
nel one at a time for ten minutes, and observed the impact 
on all other channels. For each bolometer we collected about 
660 minutes of data. 

The crosstalk coefficient between channels i and j is 
expressed as: 



Ca = AVj/AVi, 



(3) 



where Vi and Vj are the channel i and j voltages, corrected 
for thermal drift. The crosstalk matrix and a histogram of 
crosstalk levels are shown in Fig. 8. The crosstalk is mostly 
confined to nearest neighbours in the belt, channels whose 
wiring is physically close. The measured crosstalk level is 
in good agreement with ground measurements, typically 
< — 70 dB, and thus meets the requirement. A few of the po- 
larization sensitive bolometer pairs show a crosstalk around 
-60 dB. 

In the next section we see that crosstalk measurement 
from glitches shows a much lower level of crosstalk. These 



results suggest that the CPV test has measured electrical 
crosstalk in current which is unrelated to the scientific sig- 
nal. 



3.4.2. Measurements using glitches 

We used high energy glitches in one channel to study the 
impact on the signal of surrounding channels. Thousands of 
glitch events are collected for one channel, and the signals 
of all other channels for the same time period are stacked. 
The crosstalk in volts for individual glitches is defined as: 



AVj/AVi 



(4) 



where V. is the glitch amplitude in volts in the channel hit 
by a cosmic ray, and Vj the response amplitude of another 
channel j. Then, for a pair of channels i and j, the global 
voltage crosstalk coefficient is 



C ti = median(c^) 



(5) 



For SWB channels, in contrast with the CPV previous 
results, no evidence of crosstalk is seen, with an upper limit 
of — 100 dB. There are outliers in galactic channels because 
of incorrect glitch flagging. A second analysis using planet 
crossing data instead of glitches gave the same results. 

Concerning the coupling between PSB pairs, we see 
crosstalk around — 60 dB, in agreement with the CPV tests; 
however, this is likely an upper limit because it includes 
the effects of coincident cosmic ray glitches which produce 
a similar effect but are not crosstalk. 



4. Beams and time response 

4.1. Measurement of Time Response 
4.1.1. Introduction 

The time response of HFI describes the shift, in amplitude 
and phase, between the optical signal incident to each de- 
tector and the output of the readout electronics. The re- 
sponse can be approximated by a linear complex trans- 
fer function in the frequency domain. The signal band 
of HFI extends from the spin frequency of the spacecraft 
(/ S pi n ~ 16.7 mHz) to a cutoff defined by the angular size 
of the beam (14-70 Hz; see Table 4 from Lamarre et al. 
(2010)). For the channels at 100, 143, 217, and 353 GHz, 
the dipole calibration normalizes the time response at the 
spin frequency. To properly measure the sky signal at small 
scales, the time response must be characterized to high pre- 
cision across the entire signal band, spanning four decades 
from 16.7 mHz to - 100 Hz. 

The time response of bolometers typically is nearly flat 
over a signal band from zero frequency to a frequency de- 
fined by the bolometer's thermal time constant, and then 
drops sharply at higher frequencies. For the HFI bolome- 
ters, the thermal frequency is 20-50 Hz (Lamarre et al. 
2010; Holmes et al. 2008), as noted in Lamarre et al. (2010) 
and Pajot et al. (2010), however, the time response of HFI 
is not flat at very low frequencies, but exhibits a low fre- 
quency excess response (LFER). 

We define the optical beam as the instantaneous direc- 
tional response to a point source. Any sky signal is con- 
volved with this function, which is completely determined 
by the optical systems of HFI and Planck. 
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Since Planck is rotating at a nearly constant rate and 
around the same direction, the data are the convolution of 
the signal with both the beam and the time response of 
HFI. We separate the two effects and deconvolve the time 
response from the time ordered data. This deconvolution re- 
sults in a flat signal response, but necessarily amplifies any 
components of the system noise that are not rolled off by 
the bolometric response. This amplified noise is supressed 
by a low-pass filter (Planck HFI Core Team 2011b). 

4.1.2. TF10 model 

The main ingredients of the time response are: (i) heat 
propagation within the bolometer; (ii) signal modulation at 
a frequency of / mo d = 90.188 Hz performed by reversing the 
bolometer bias current; (iii) the effect of parasitic capaci- 
tance along the high impedance wiring between the bolome- 
ter and the first electronics stage (JFETs); (iv) band-pass 
filtering, to reject the low frequency and high frequency 
white noise in the electronics; (v) signal averaging and sam- 
pling; and (vi) demodulation. 

Because of the complexity of this sequence, a phe- 
nomcnological approach was chosen to build the time re- 
sponse model. The time response is written as the product 
of three factors: 



Hia(f) — -ffboio x H ICS x H, 



filter 



(6) 



Schematically, the first factor takes into account step 
(i), the second factor describes a resonance effect that re- 
sults from the combination of steps (ii) and (iii) , while the 
purpose of -f/fiitcr is to account for step (iv). 

Detailed analysis and measurements of heat propaga- 
tion within the bolometer have shown that -Hboio is given 
by the algebraic sum of three single pole low pass filters. 
Explicitly: 



bolo 



»=1,3 



with 6 parameters (oi, a^, (13, ti, T2, T3) 

1+ P7 (2tt/) 2 
rcs l-p s (27rfy+jp g (2nf) 

with 3 free parameters (j»7,pg,pg), 

1 



filter 



(f/F mod ) 2 



l-p 10 (27T/) 2 



J (//^filter) 2 



(7) 



(8) 



(9) 



with one free parameter (pro)- A total of 10 free parameters 
describe this model, as indicated by its name. See Fig. 9 for 
an illustration of the three components of the time response 
model TF10 for a typical 217 GHz channel. 

The parameter Ffiitcr characterizes the rejection filter 
width and is kept fixed to 6 Hz in the fitting process. Besides 
the fact that this phcnomcnological model is physically mo- 
tivated, this parameterization: 

— ensures causality 

— satisfies H(-f) = H*(f) 

— goes to 1 when / goes to zero (because we define a\ + 
a-2 + 0-3 = 1), while it goes to when / goes to infinity 

— includes enough parameters to provide the necessary 
flexibility to fit the time response data of all 52 bolome- 
ters. 




40 60 
Frequency(Hz) 

Figure 9. The amplitude and phase (in radians) of the 
three components of the TF10 model of the time response. 
The solid blue line is i?boio(/), the dotted green line shows 
-ff filter (/) an d the dashed red line shows H res . The solid 
black line is Hiq(J), the product of the three components. 
The vertical dotted black line shows the signal frequency 
where the beam of a 217 GHz channel cuts the signal power 
by half. 



4.1.3. Fitting the TF10 Model to Ground Data 

To obtain the 10 x 52 parameter values, we used three sets 
of pre-launch measurements, (i) The bolometer response 
was measured at 10 different frequencies by illuminating all 
52 bolometers with a chopped light source, (ii) Other mea- 
surements were done using carbon fibers as light sources; 
the latter were alternately turned on and off at a variable 
frequency, (iii) The bolometer bias currents were periodi- 
cally stepped up and lowered by a small amount. By adding 
a square wave to the DC current, temperature steps are in- 
duced, simulating turning on and off a light source (the 
analysis of these data requires bolometer modelling) . 

None of these measurements was absolutely normalized; 
all compared the relative response to inputs of various fre- 
quencies. While measurement (i) only provided the ampli- 
tude of the time response, measurements (ii) and (iii) pro- 
vided both the amplitude and the phase. Note that for the 
phase analysis, because of the lack of precise knowledge of 
the time origin t of the light/current pulses, a fourth fac- 
tor exp(j27rAio/) is introduced in the expression of Hiq(J), 
where the additional parameter Atg represents the uncer- 
tainty in time. 

Among the three sets of measurements, the carbon fiber 
set covered the largest frequency range. Thus it was the best 
for building the transfer function model described above 
and for investigating its main features. However, it involves 
uncertainties that could be resolved only with a very de- 
tailed simulation of the set-up, therefore it was not used to 
calculate the final set of parameter values. The final values 
were calculated from data sets (i) and (iii) , whose frequency 
ranges are complementary, 2-140 Hz and 0.0167-10 Hz, re- 
spectively. Since no absolute normalization was available, 
the two datasets were matched in the overlap frequency 
range. 
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The fitting of the analytic expression given above to 
the merged data was done in the range between 16.7 mHz 
and 120 Hz. The 52 fits have a % 2 /DoF distribution whose 
mean value is 1.13, indicating that the model is adequate 
to describe the data. The numerical values of pio displayed 
a small spread, c moan < 6 x 10~ 4 . This parameter was set 
at its mean value to calculate the 52 covariance matrices of 
the nine remaining parameters, which were useful in prop- 
agating the statistical errors. 

As described below, the time response thus obtained 
was further tuned and checked using in-flight observations, 
in particular signals produced by planets and by cosmic 
rays (glitches). 

An alternative model has also been defined, based on the 
analytical expression of the steps (ii) to (vi) of section 4.1.2. 
Based on a closer analysis of the electronics stages, this 
model is more physically motivated than the 10-parameter 
model. It requires only 8 parameters and provides better 
results near the modulation frequency. Nevertheless the 
model has not been used in the current data release. It 
is only used as a benchmark, to check possible systematic 
effects in the current release. Most of the effects of the dif- 
ference between the models disappear when the data are 
low-pass filtered. 

4.1.4. Fitting TF10 to Flight Data 

The planets Mars, Jupiter, and Saturn are bright, compact 
sources that are suitable for measuring the beam and pro- 
vide a near-delta-function stimulus to the system that can 
be used to constrain the time response. During the first sky 
survey, Planck observed Mars twice and Jupiter and Saturn 
once (Planck HFI Core Team 2011b). During a planet ob- 
servation, the spacecraft scans in its usual observing mode 
(Planck Collaboration 2011a), shifting the spin axis in 2' 
steps along a cycloidal path on the sky. Since planets are 
close to the ecliptic plane, the coverage in the cross-scan di- 
rection is not as fine as in the scan direction. In the case of 
Jupiter and Saturn, each channel observes the planet once 
per rotation for a period of approximately 6 hours (9 peri- 
ods of stationary pointing, or "rings"). Because Mars has 
a large proper motion, the first observation lasted 12 hours 
(or 18 rings). 

We use the forward-sense time domain approach 
(Huffenberger et al. 2010) to simultaneously fit for Gaussian 
beam parameters and TF10 time response parameters. A 
custom processing pipeline avoids filtering the data. We 
extract the raw bolometer signal and demodulate it using 
the parity bit. We use the flags created by the time ordered 
information (TOI) processing pipeline to exclude data sam- 
ples contaminated by cosmic rays, and we additionally flag 
all data samples where the nonlinear gain correction is more 
than 0.1%. We use Horizons 2 emphemerides to compute the 
pointing of each horn relative to the planet center. 

The time domain signal from the planet is modeled as an 
elliptical Gaussian convolved with the TF10 time response 
as follows: 



d(t) = H 10 * A(t)G [x(t); x , e, 6 FWHM , i>] 



(10) 



where the Gaussian optical beam model G is parameterized 
as in Eqs. 9-11 of Huffenberger et al. (2010), except the 
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Figure 10. Comparison of the impulse response of chan- 
nel 143-2a (red curve) and a template made from stack- 
ing glitch events (black curve). Noise begins to dominate 
further in the timeline. The ringing observed at the mod- 
ulation frequency is generated by the electronics rejection 
filter. 



planet amplitude is parameterized with a disk temperature 
rather than a single amplitude: 



A(t) = T disk 



n p {t) 



(ii) 



ssd. jpl .nasa. gov/?horizons 



where Tdisk is the whole-disk temperature of the planet, 
rip is the solid angle of the planet, which can vary signifi- 
cantly during the observation, and flj, is the solid angle of 
the beam. Q p is computed using Horizons, which is pro- 
grammed with Planck's orbit. 

The free parameters of the fit are the six parameters 
of the time response corresponding to -ffboio, the two com- 
ponents of the centroid of the beam Xq, the mean FWHM 
$fwhm, the cllipticity e, the ellipse orientation angle ip, and 
the planetary disk temperature Tdisk- The four parameters 
describing the electronics are somewhat degenerate with 
the bolometer part of the time response, and we fix them 
at the ground-based values. 

Because of the large nonlinear response and highly non- 
Gaussian beams at 545 and 857 GHz, we do not perform 
fits to the planet data at these frequencies. Instead we rely 
on pre-launch fits for the time response. 

By taking the Fourier transform of the time response 
function derived on planets, one obtains the system re- 
sponse to a Dirac impulse. This response can be compared 
to the glitches generated by cosmic rays that deposit energy 
in the sensor grids. 

The glitches detected by HFI are sampled with time 
steps l/(2F mo( i). However, the glitches can be superresolved 
in time by normalizing, phasing, and stacking single glitch 
events (Crill et al. 2003). This gives glitch templates for 
each channel (Alexandre Sauve, private communication) 
that are effectively sampled at a much higher frequency. 

Figure 10 shows the comparison between a superre- 
solved glitch template and the corresponding calculated re- 
sponse. There is good agreement in general, but there are 
discrepancies at high frequency (/ > 100 Hz). The physical 
model for the electronics transfer functions briefly described 
at the end of section 4.1.3 suppresses this discrepancy at 
high frequency. 
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Each planet observation suffers unique systematic ef- 
fects, so a comparison of the time response recovered on 
each gives a good assessment of the effects of these different 
systematics. Mars has a large proper motion, giving excel- 
lent sampling in the cross-scan direction; however, there is 
a known diurnal variability in its brightness temperature 
(Swinyard et al. 2010). Jupiter has a large angular diam- 
eter (48") relative to the HFI beam size, and Saturn and 
Jupiter are so bright that the HFI detectors are driven sig- 
nificantly nonlinear (see Table 3). Nonetheless, we find that 
the time response is consistent to 0.1-0.5% when recovered 
from each of the planets individually as well as from all 
planets simultaneously. 

A further cross-check is done by stacking planet scans 
to build a superresolved planet timeline. Time response pa- 
rameters are fit to the superresolved planet using the as- 
sumption of a near-Gaussian beam profile, and are consis- 
tent with the first approach. 

The in-flight time response differs from the ground- 
based time-response by at worst 1.5% between 1Hz and 
40 Hz. We do not include this difference in the final er- 
ror budget, because it is likely that the time response has 
changed due to differences in background conditions. 

4.1.5. Low Frequency Excess Response 

The HFI bolometers show low frequency excess response 
(LFER), (Lamarre et al. 2010). Though the planets are 
bright, the short impulse they provide is close to a delta 
function and the energy is spread evenly across nearly all 
harmonic components. In combination with low frequency 
noise, the measurements are not sensitive to frequencies 
below ~ 0.5 Hz; so with planet observations alone we can- 
not constrain an excess response at very low frequency. We 
maintain the ground measurements as our best estimation 
of the LFER. In the ground-based measurements the bias 
step vs. carbon fiber differs by at most 1.5% at low fre- 
quency, so we assign a systematic error of 1.5% for frequen- 
cies below 0.5 Hz. 

For future data releases, we will use the difference of sky 
signal between surveys to constrain the LFER. 

4.1.6. Summary of Errors in the Time Response 

As noted above, the data represent the combined effect of 
the time response and the optical beam. The time response, 
however, is not degenerate with a Gaussian parameteriza- 
tion of the beam; the true beams deviate from a Gaussian 
shape at the several percent level near the main lobe, while 
time response effects tend to give the beam an extended tail 
following the planet in the scan direction. The Gaussian as- 
sumption could slightly bias the recovered time response; 
however, any residual bias is captured in the measurement 
of the post-deconvolution scanning beam (Planck HFI Core 
Team 2011b). 

Because of the high signal-to-noise ratio of the planet 
data, statistical errors in the fit are small, so we assess the 
systematic errors in the resulting time response by checking 
the consistency of various methods of recovering the time 
response. We fit to different combinations of planet data: 
Mars, Jupiter, and Saturn data separately and all of the 
data simultaneously to check for systematics resulting from 
various planets. Additionally we compare the planet-fitted 



time response with ground-based data and with the impulse 
response from cosmic ray glitches. 

Our final error budget is as follows: 

— Low frequency (/ < 0.5 Hz): the errors are dominated 
by the possibility of a low frequency excess response 
below 0.5 Hz at a level of 1.5%. 

— Middle frequency (0.5 Hz < / < 50 Hz): We set an error 
bar between 0.1% and 0.5% depending on the channel. 
This error bar is set by the consistency in results from 
different sets of planet data. 

— High frequency: (/ > 50 Hz) Our empirical model of 
the electronics in the TF10 model does not describe 
the system very well at these frequencies, as shown by 
some disagreement between the glitches and the TF10 
impulse response. However, for this data release, the 
low-pass filter applied to the data and the beam cutoff 
reduce the importance of this frequency band. 

The Planck scan strategy is such that the same region of 
the sky is observed scanning in nearly opposite directions 
six months apart. An error in the time response is high- 
lighted in the difference of maps obtained from the first 
six months and the second six months of the survey. This 
difference map shows some level of contamination, in par- 
ticular near the Galactic plane, where the signal is higher. 
The same level of contamination is observed in simulations 
in which the data are generated with a transfer function, 
and analysed with a different one, in order to mimic the 
uncertainties described above. With this technique, we val- 
idate the error budget. 

4.2. Optical Beams 

The optical beam is defined (Sect. 4.1.1) as the instanta- 
neous directional response to a point source. For HFI, the 
optical beams for each channel are determined by the tele- 
scope, the horn antennas in the focal plane and, for the po- 
larized channels, by the orientations of their respective po- 
larization sensitive bolometers (PSBs) (Maffei et al. 2010). 
Model calculations of the beams are essential, since it was 
possible to measure only a limited number of beams in 
the telescope far-field before launch. The 545 and 857 GHz 
channels, which employ multimoded corrugated horns and 
waveguides, were not included in this campaign. (The opti- 
cal beam is related to, but is not the same as the scanning 
beam defined in Planck HFI Core Team (2011b) and used 
for data analysis purposes.) Tauber et al. (2010b) reported 
the best pre-launch expectations for the optical beams, ob- 
tained using physical optics calculations with CORRUG 3 
and GRASP 4 . Table 4 compares the calculated and mea- 
sured (Sect. 4.1.2) beams for the single-moded channels (up 
to 353 GHz). 

For these channels the pre-launch calculations of 
FWHM and cllipticity and measured Mars values agree to 
within a few percent. These differences are contained within 
2.7<t of the data errors. The main source of discrepancy 
could be a slight misalignment of the pre-launch telescope 
model with respect to the actual in-flight telescope geom- 
etry, which is currently being investigated (Jensen & et al 
2010). 

3 SMT Consultancies Ltd. www.smtconsultancies.co.uk 
TICRA, www.ticra.com 
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Table 4. Comparison of pre-launch calculations and measured parameters for the HFI optical beams (band averages) . 
Standard deviations a are computed as the dispersion between the Saturn, Jupiter, and Mars data for each given channel. 
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Figure 11. The distribution of the HFI beams on the sky 
relative to the telescope boresight as viewed from infinity. 
Contours of the Gauss-Hermite decomposition of the Mars 
data at 1%, 10%, and 50% power levels from the peak. 
For the photometers containing a pair of PSBs, the average 
beam of the two PSBs is shown. 



Table 5 reports our best knowledge of the FWHM of the 
optical beams for each channel. We stress that this table 
does not provide parameters of the scanning beam of the 
processed data, which accounts for the additional effects 
of the instrument time response and of the time domain 
filtering in the data processing (Planck HFI Core Team 
2011b). 

As reported in Maffei et al. (2010), the 545 GHz and 
857 GHz channels are multimoded (more than one elec- 
tromagnetic mode propagating through the horn anten- 
nas) and their optical beams are markedly non-Gaussian. 
The understanding of these channels through simulations 
has progressed since Planck was launched, especially in 
the characterization of their modal content (Murphy et al. 
2010). 

In Table 4 we compare pre-launch calculations of the 
beams with the beams measured with Mars. Differences in 
FWHM are less than 7%. We stress that this discrepancy 
does not impact the scientific products of the Planck mis- 






Figure 12. The "dimpling effect" as seen at 545 GHz (left 
panel) and 857 GHz (right panel). The grid spacing is 10'. 
The color scale is in dB normalized to the peak signal of 
Jupiter. 



sion since the scanning beams are the ones to be used for 
data analysis purposes. From an instrumental point of view, 
the inflight measurements must obviously be considered as 
the reference for the performance of these channels. 

The development of the HFI multimoded channels ne- 
cessitated the novel extension of previously existing mod- 
elling techniques for the analysis of the corrugated horn 
antennas and waveguides, as well as for the propagation 
of partially coherent fields (modes) through the telescope 
onto the sky (Murphy et al. 2001). Extensive pre-launch 
measurement campaigns were conducted for all the HFI 
horn antenna/filter assemblies (Ade et al. 2010). The HFI 
multimoded channels are suitable for the scientific goals of 
Planck. Nevertheless, for future instruments, more research 
can be envisaged in this field. The characterization of the 
modal filtering in the horn-waveguide assembly and the un- 
derstanding of the coupling of the waveguide modes to the 
detector need further theoretical and experimental study. 

The similarity of the pre-launch expectations to our cur- 
rent knowledge of the HFI focal plane (beams and their 
positions on the sky) tells us that the overall structural in- 
tegrity of the focal plane has been preserved after launch. 
Furthermore, the optical beams as measured on Mars are 
shown in Fig. 11 and can be compared with the equivalent 
representations of the focal plane layout based on calcu- 
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lations in earlier papers (Maffei et al. 2010; Tauber et al. 
2010b). A detailed account of the full focal plane recon- 
struction can be found in Planck HFI Core Team (2011b). 

There is a "dimpling" of the reflector surfaces from the 
irregular print-through of the honeycomb support struc- 
tures on the reflector surfaces themselves (Tauber et al. 
2010b). GRASP calculations predict that this will gener- 
ate a series of rings of narrow bright grating lobes around 
the main lobe. Since the small-scale details of the dimpling 



structure of the Planck reflectors are irregular, these grat- 
ing lobes tend to merge with the overall power scattered by 
the reflector surfaces (Ruze scattering (Ruze 1966)). Fig. 12 
shows a HEALPIX (Gorski et al. 2005) map of the first sur- 
vey observation of Jupiter minus the second survey observa- 
tion of the same sky region to remove the sky background. 
We see the first ring of grating lobes as expected in the 
map from all 545 and 857 GHz channels, where the signal- 
to-noise ratio on the planets is highest. The inner 15' of 
the beam is saturated and does not appear in the map. At 
857 GHz, the discrete grating lobes appear at level below 
— 35 dB with respect to the peak (~ 30 dB), and represent a 
negligible fraction of the total beam throughput. The shoul- 
der of the beam, extending radially to ~ 15', represents a 
larger contribution to the throughput, ranging from 0.5% 
to a few percent for the CMB and sub-mm channels, re- 
spectively. 

5. Noise properties 

The Planck HFI is the first example of space-based bolome- 
ters, continuously cooled to 100 mK for several years. 
Although the detectors were thoroughly tested on the 
ground (Lamarre et al. 2010; Pajot et al. 2010), it remained 
then to be seen how they would behave in the L2 space envi- 
ronment. We describe here the noise properties of the HFI 
in the first year of operation, focusing on the differences 
between space and ground performance. 

This section deals with the Gaussian part of the noise. 
Section (6) describes the systematic effects that have been 
analyzed in the data so far. 

An example of raw time ordered information (TOI) is 
shown in Fig. 13. The TOI is dominated by the signal 
from the CMB dipole, Galactic emission, point sources, and 
glitches. Therefore, the noise properties cannot be directly 
deduced from the TOI. We first describe the general method 
used to evaluate the noise, then we give general statements 
on the noise properties. 

5.1. Noise estimation 

The Level-2 detnoise pipeline (Planck HFI Core Team 
2011b) is used to determine a noise power spectrum, from 
which one extracts the noise equivalent power (NEP) of 
the detectors (see Planck HFI Core Team (2011b) for a full 
description). The pipeline uses redundancies in the obser- 
vations to determine an estimate of the sky signal, which is 
then subtracted from the full TOI to produce a pure noise 
timeline. The signal estimates are the integration of typ- 
ically 40 circles of data at a constant spin axis pointing. 
The average signal, binned in spin phase, provides an accu- 
rate estimate of the signal. This signal as a function of spin 
phase is then subtracted from the TOI. The residual is an 
estimate of the instantaneous noise. Power spectra of this 
residual timeline are then obtained for each pointing period 
(see Fig. 14) and fit for the white noise level, i.e., the NEP, 
in the spectral region between 0.6 and 2.5 Hz. The lower 
limit of 0.6 Hz is high enough that the low frequency excess 
noise can be neglected and the upper limit small enough to 
keep the time response near to its value at low frequencies 
(16mHz) at which the instrument is calibrated. 

The noise is stable at a level better than 10% in the 
majority of detectors. Exceptions are: (1) a few rings 
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Figure 13. Examples of raw (unprocessed) TOI for one bolometer at each of six HFI frequencies and one dark bolometer. 
Slightly more than two scan circles are shown. The TOI is dominated by the CMB dipole, the Galactic dust emission, 
point sources, and glitches. The relative part of glitches is over represented on these plots due to the thickness of the 
lines that is larger than the real glitch duration. 



with unusual events that contaminate the measurement, 
e.g., poorly corrected/flagged glitches or passage over very 
strong sources such as the Galactic centre, especially at 
857 GHz; (2) a weak trend, smaller than 1% in amplitude, 
that correlates with the duration of the pointing period 
(an expected bias due to the ring average signal removal); 
(3) bolometers affected by random telegraph signals (RTS) 
(see next Section); and (4) some uncorrelated jumps in the 
noise levels for about ten bolometers at the 30% level for 
isolated periods of a few days. The overall result is that a 
very clear baseline value can be identified and can be used 
to determine the NEP of each bolometer. This is then con- 
verted to NEAT with the help of the flux calibration. The 
NEATs thus measured are given in Table 6. The quoted 
uncertainties are derived from the rms of the NEPs in a 
band around the baseline. 



noise filtered by the time response of the bolometer, 
the readout electronics, and the TOI processing. 

— Electronics and Johnson noise, which produce noise that 
is nearly white across the frequency band, but with a 
sharp decrease at the high frequency end due to the 
on-board data handling and the TOI filtering. 

— The 4K lines (Sect 6), appearing as residuals in the 
spectra. 

— The energy deposited by cosmic rays on the bolome- 
ters, which appears as "glitches", i.e., positive peaks 
in the signal, which are removed by the TOI processing 
(Sect. 6, and Planck HFI Core Team (2011b)). Residuals 
from glitches appear in the noise spectrum as a bump 
between 0.1 and 1Hz. 

— Low frequency excess (LFE) noise, which is present be- 
low about lOOmHz. 



. The last three sources of noise are detailed in Section 6. 

5.2. The noise components „, , , r , , , r „ r0/ , 

r ihere is additional noise (of the order ol 0.5% or less) 

The detector noise is described by the combination of sev- due to the on-board quantization of the data before trans- 

eral components: mission. In general, the noise level, as measured by the 

NEP, is between 10 and 20aWHz~ 1 / 2 for the 100 to 

— Photon and bolometer noise, which appear as white 353 GHz channels, and between 20 and 40aWHz -1 / 2 for 
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Figure 14. Typical power spectrum amplitude of bolometers 143-5 and 545-2. For the upper panels, this is the power 
spectrum density of valid samples, after an average ring (the sky signal) has been subtracted from the TOI. Stacking of 
the result for 200 rings is shown in the lower panel. Here, the instrument time response is not deconvolved from the data. 



the 545 and 857 GHz channels. It is in line with the ground- 
based expectations and the lower estimate of the back- 
ground load with a detector-to-detector variability of less 
than 20% (see Sect. 8). 

Due to the AC bias modulation scheme, the 1// noise 
from the electronics is aliased near the modulation fre- 
quency where it is heavily filtered out. The benefit of this 
scheme is visible on the noise power spectrum of the 10 Mil 
resistor which shows a flat spectrum at the Johnson value 
down to 1 mHz, a tribute to the electronic chain stability. 

At the present time, we assume that the LFE noise, not 
observed in ground-based measurements, is mostly due to 
the 100 mK bolometer plate fluctuations. While drifts in 
the 100 mK stage that are correlated between bolometers 
are removed, there are likely local temperature fluctuations 
due to particle energy deposited close to each detector. 

6. First assessment of systematic effects 

6.1. 4 K lines 

The fundamental frequency of the 4 He-JT cooler 
(40.083373 Hz) is phase-locked to the frequency of the data 
acquisition (180.37518 Hz) in a ratio of 2 to 9. EMI/EMC 
impacts the TOI only as very narrow lines. Unfortunately, 
in flight, unlike in ground-based measurements, these lines 
are not stable. The 4K line variations are illustrated in 



Fig. 15. The variability of the lines is in part due to tem- 
perature fluctuations in the service module of the Planck 
spacecraft. Indeed, some of the variability was related to the 
power cycling of the data transponder which, for stability 
reasons, has been kept on continuously since 25 January 
2010 (OD 258, Planck Collaboration 2011b, see Fig. 16). 



6.2. Abnormal noise in the electronics 



Of the 54 bolometers on HFI, three show a significant RTS, 
also known as "popcorn noise." These are 143-8, 545-3, and 
857-4. Fig 17 illustrates their behaviour. The noise timeline 
clearly exhibits a two-level system. The three RTS bolome- 
ters in flight are the ones where RTS occurred most fre- 
quently in ground measurements. However, in flight: 1) the 
level difference is well above the noise (at least ten times the 
rms); 2) the two-level system can be a three-level system 
or even larger; and 3) the RTS is intermittent. For large 
duration, it can be unnoticeable, especially for the 857-4 
bolometer. 

In an unrelated fashion, we see uncorrelated jumps in 
the noise TOI of many bolometers at a rate of just a few 
every year. 
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Figure 16. Zoom on the four main 4K line systematics amplitude measured in the TOI processing on the test resistance. 
{Left) Period when the transponder was switched on once per day for the 3-hours Daily Tele-Communication Period 
(DTCP). (Right) Period when the transponder was kept on at all times, later in the mission. 




Figure 15. Typical trend of the cosine and sine coefficient 
variation of the four main 4 K lines measured in the TOI 
processing on the test resistance at 10, 30, 50, and 70 Hz 



6.3. Cosmic rays and their effects 

Energy is deposited by cosmic rays in various parts of the 
HFI instrument. We observe these events in the TOIs of 
all detectors as a signal peak characterised by a very short 




Figure 17. Random telegraphic signal in the noise time- 
line (fW) of bolometer 545-3 plotted vs. time. RTS is here 
a two-level signal with random occurrences of flipping. Full 
sampling is in black dots. A smoothed version (by 41 sam- 
ples) is plotted with a red line. 



rise time (less than 1.5 ms) and an exponential decay. These 
events are called glitches. The other effect of the cosmic rays 
is a thermal input to the bolometer plate, which induces 
low frequency noise on the bolometers. Thermal effects are 
described in Planck HFI Core Team (2011b) and their very 
long term consequences detailed in Sect. 7. 
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Figure 18. SREM hit count and HFI bolometer average 
glitch rate evolution. SREM TCI hit counts measure the 
protons with a deposited energy larger than 0.085 MeV. 

6.4. Cosmic ray-induced glitch spectrum seen by Planck 

Cosmic rays consist of two main components at the L2 lo- 
cation of Planck: the Solar component and the Galactic 
component. The Solar component is at low energy (a few 
keV) except during flares, when energies can reach GeV. 
HFI is immune to the low energy component and no ma- 
jor flares have yet been recorded. The Galactic component 
(Shikaze et al. 2007) with a maximum between roughly 
300 MeV and 1 GeV, is modulated by Solar activity. The 
Planck mission began during the weakest solar activity for 
a century (McDonald et al. 2010). Hence, the Galactic cos- 
mic ray component is expected to be at its highest level. 

The glitch rate evolution (Fig. 18) follows closely the 
proton monitoring by the space radiation environnement 
monitor (SREM), a piggy-back experiment mounted on the 
Planck spacecraft. This figure shows that cosmic rays are 
the main source of HFI glitches. The glitch rate tended to 
decrease during the mission since January 2010 due to the 
slow increase of the Solar activity. 

This glitch rate can be understood as the sum of two 
interaction modes, depending on whether the cosmic ray 
has a direct or indirect interaction with the bolometer. High 
energy cosmic rays can also interact with the bolometer 
plate and induce thermal effects and correlated glitches on 
the bolometers. These are dealt with in Sect. 7. The glitch 
characteristics also depend on the location of the energy 
deposit within the bolometer: the thermistor, the absorbing 
grid, or the bolometer housing. 

6.4.1. Direct interaction 

Cosmic ray particles can deposit energy directly on the 
thermistor or the absorbing grid. This is observed at well- 
defined deposited energies corresponding to the thickness 
of the element. Typically thermistor hits are about 20keV, 
whereas grid hits are about 2 keV. These events occur at a 
rate of a few per minute. 

6.4.2. Indirect interaction 

Cosmic ray particles can also deposit energy indirectly. All 
particles crossing some matter produce a shower of sec- 
ondary electrons, through ionization, that are mostly ab- 
sorbed in the matter nearby. However, interactions occur- 
ring within microns of the internal surface of the bolometer 
box produce a shower of free secondary particles. A frac- 
tion of these particles is absorbed by the thermistor and the 
grid of the bolometer. This explains the large coincidence 



rate of gliches between PSB bolometers a and b sharing 
the same mounting structure. The energy of those glitches 
follows a power law distribution spanning the whole range, 
from the detection threshold to the saturation level. This 
spectrum is expected for the delta and secondary electrons 
produced via the ionization process. The total rate of these 
events is typically one per second, and thus dominates the 
total counts shown in Fig. 19. 

A more detailed description of the effect of cosmic 
rays on HFI detectors is postponed to a dedicated paper, 
and glitch handling in the data processing is described by 
Planck HFI Core Team (2011b). 

7. Instrument stability 

The radiative power reaching each bolometer is the co- 
addition of the flux from the sky and of the thermal emis- 
sion of all optical elements "seen" from the detector: filters, 
horns, telescope reflectors, shields, and mechanical parts 
visible in the side-lobes. In addition, fluctuations of the 
heat sink temperature (the bolometer plate) appear like an 
optical signal. Any change in any of the parameters (tem- 
perature, emissivity, geometrical coefficient) driving these 
sources may be visible in the bolometer signal as a "DC 
level", i.e. a stable or very slowly varying (days) compo- 
nent in the signal. 

Monitoring the "DC level" supposes that one is able to 
remove the varying sky signal from the stable sources. This 
is done in the map-making process (Planck HFI Core Team 
2011b) by using the redundancy of the scanning strategy. 
Fig. 20 shows for the 217 GHz bolometers the history of 
the DC level during nearly one year. All follow a pattern 
similar to that of the cosmic ray activity measured by the 
SREM (see Sect. 6 and Fig. 18), which indicates that cos- 
mic rays are at the origin of the measured signal. One can 
check on families of bolometers with non-uniform heat leaks 
G that this signal is directly related to temperature vari- 
ations of the bolometer plate and not to external optical 
sources. In fact, we can see here residual fluctuations that 
the PID of the bolometer plate fails to compensate because 
its efficiency is far from one. The similarity of Fig. 20 and 
Fig. 18, also shows that the effect of gain variations and 
of DC level drifts of the readout electronics is small with 
respect to other sources of signal drifts. 

It should be noted that the "DC level" variation of 
217 GHz bolometers is equivalent to an optical power of 
a couple of femtowatts, while the total background power 
on these bolometers is about 1 pW. This fluctuation is 
mainly due to the energy deposited by cosmic rays on the 
bolometer plate, which means that the "equivalent power" 
of the other sources of temperature fluctuation and of op- 
tical background fluctuations are no more than a fraction 
of femtowatt, i.e. less than one part per thousand of the 
background. 

The change of gain induced by the DC level variations 
can be estimated from the non-linearity measurements (see 
Sect. 3.3.2). In the case considered in Fig. 20, the relative 
gain change is of the order of a few 10~ 4 . 

During the CPV phase, the readout electronics was 
"balanced", i.e. the offset parameter was tuned to get a sig- 
nal near to zero. During the first year of operation, and for 
all bolometers, deviations from this point remained small 
with respect to the total range of acceptable values. In 
consequence, no re-tuning of the readout electronics was 
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Figure 19. Glitch rate of all HFI bolometers. An average over the first sky survey has been performed. The asymmetry 
between PSB bolometers sharing the same horn is an effect of detection threshold and asymmetric time constant properties 
between PSB a and b. 
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Figure 20. The drift in the DC level of the 217 GHz SWB 

bolometers, in femto Watts of equivalent power in the de- 
tector, for the first year of operation. 



needed during this period and it is expected to be the same 
up to the end of the mission. 



8. Main performance parameters 

The primary difference between the in-flight and pre-launch 
performance of the HFI derives from the relatively high rate 
of cosmic rays in the L2 environment. At the energies of in- 
terest, the low level of solar activity results in an elevated 
cosmic ray flux. The glitches that result from cosmic ray 
events must be identified and removed from the time or- 
dered information prior to processing the data into maps. 
The TOI processing also removes a significant fraction of 
the common mode component that appears in the bolome- 
ter TOIs at low frequencies. A residual low-frequency com- 
ponent is removed during the map-making process (Planck 
HFI Core Team 2011b). 



Table 6 summarizes the noise properties of the processed 
TOI (Planck HFI Core Team 2011b), by the following pa- 
rameters: 

— A white noise model. NEPi is the average of the Noise 
Equivalent Power spectrum in the 0.6-2.5 Hz range. 

— A model with a white noise NEP2 plus a low frequency 
component: NEP = NEP 2 [1 + (f k nee/f) Q ] 

— The sensitivity NEAT^mb to temperature differences 
of the CMB . Note that this quantity is not particularly 
relevant for the channels at 545 and 857 GHz for which 
it takes large values that are highly dependent on the 
details of the spectral transmission for each detector. 

— The sensitivity NEAT rj to temperature differences for 
sources observed in the Rayleigh-Jeans regime. 

Figure 21 compares the goal, pre-launch and in-flight 
NEATs. The average in-flight NEATs are 27% higher than 
the pre-launch NEATs. While the pre-launch and in-flight 
values are not directly comparable due to differences in the 
processing, these differences can account for less than half 
of the observed variation. The remaining part is attributed 
to residual contamination from cosmic rays that are not 
completely removed in the current TOI processing. 

The sensitivity goals are taken from Table 1 of Lamarre 
et al. (2010), and are consistent with Table 1.3 of Planck 
Collaboration (2005) corrected for the use of PSBs at 
100 GHz. Note that Fig. 21 supersedes Figure 11 of Lamarre 
et al. (2010) in which requirements and goals were im- 
properly plotted. The in-flight sensitivities estimated from 
NEPi exceed the goals, which are defined by a total noise 
level equal to twice the expected contribution of photon 
noise. The average measured NEPi is typically 70% of the 
initial goal. The improvement in the NEP over the design 
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Figure 21. Noise Equivalent Delta Temperature measured 
on the ground and in-flight with slightly different tools 

goals is primarily the result of having reduced the photon 
background through careful design of the optical system. 



plate that induce a common mode component to the sci- 
entific signal and the focal plane thermometry Although 
a component correlated with the dark bolometer outputs 
is removed during the TOI processing, a residual low fre- 
quency contribution is removed at the map-making stage. 
With the exception of the three detectors affected by tele- 
graph noise, the sensitivity measured above 0.6 Hz exceeds 
the design goals of all HFI channels. After the removal of 
the residual low frequency noise, the final sensitivity of the 
frequency maps exceeds the mission requirements, and ap- 
proaches the goals, as described in the companion paper 
Planck HFI Core Team (2011b). 
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indicating which technical or scientific activities they have been in- 
volved in, can be found at www . rssd. esa. int/index . php?project= 
PLANCK&page=Planck_Collaboration. 



9. Conclusions 

We report on the in-flight performance of the High 
Frequency Instrument on board the Planck satellite. These 
results are derived from the data obtained during a ded- 
icated period of diagnostic testing prior to the initiation 
of the scientific survey, as well as an analysis of the sur- 
vey data that form the basis of the early release scientific 
products. 

With the exception of a single anomaly in the operation 
of the 4 He-JT cooler, the HFI has operated nominally since 
launch. The settings of the readout electronics determined 
during pre-launch testing were found to be very near the 
optimal value in flight and were applied without any mod- 
ification. A random telegraphic signal is observed in the 
same three channels that exhibited this behaviour during 
the final pre-launch testing. These channels are currently 
excluded from the scientific analysis. The instrument op- 
eration has been extremely stable during the first year of 
operation, requiring no adjustment of the readout electron- 
ics. 

The optical design, and the alignment of the optical as- 
sembly, relied on both theoretical analysis and testing at 
the subsystem level. The beams of the 545 and 857 GHz 
channels, which employ multimoded corrugated horns and 
waveguide, could not be measured on the ground. The ac- 
tual beam widths of these channels measured on planets 
are in general smaller than the design goals and estimated 
values. The optical properties of the single mode channels 
are in excellent agreement with the design expectations. 

A higher than expected cosmic ray flux, related to the 
level of Solar activity, results in a manageable loss of sig- 
nal and degradation of thermal stability. Discrete cosmic 
ray events result in glitches in the scientific signal that 
are flagged and removed by an algorithm making use of 
the signal redundancy in the timeline. In addition to these 
single events, the cosmic ray flux contributes a significant 
thermal load on the sub-kelvin stage. Variations in this 
flux produce low-frequency fluctuations in the bolometer 



References 

Ade, P. A. R., Savini, G., Sudiwala, R., et al. 2010, A&A, 520, A11+ 
Benoit, A., Ade, P., Amblard, A., et al. 2003a, A&A, 399, L25 
Benolt, A., Ade, P., Amblard, A., et al. 2003b, A&A, 399, L19 
Benoit, A., Sirbi, A., Bradshaw, T., et al. 1997, in ESA Special 
Publication, Vol. 400, Sixth European Symposium on Space 
Environmental Control Systems, ed. T.-D. Guyenne, 497 — f 
Bersanelli, M., Mandolesi, N., Butler, R. C, et al. 2010, A&A, 520, 
A4+ 

Bhandari, P., Bowman, R. C, Chave, R. G., et al. 2000, Astrophysical 

Letters Communications, 37, 227 
Bhandari, P., Prina, M., Bowman, R. C, et al. 2004, Cryogenics, 44, 

395 

Bock, J. J., Chen, D., Mauskopf, P. D., & Lange, A. E. 1995, Space 

Science Reviews, 74, 229 
Bradshaw, T. W. & Orlowska, A. H. 1997, in ESA Special Publication, 

Vol. 400, Sixth European Symposium on Space Environmental 

Control Systems, ed. T.-D. Guyenne, 465 — f 
Catalano, A. 2008, PhD thesis, Ecole doctorale Astronomie et 

Astrophysique d'lle de France 
Catalano, A., Coulais, A., & Lamarre, J. 2010, Appl. Opt., 49, 5938 
Church, S. E., Philhour, B. J., Lange, A. E., et al. 1996, in ESA 

Special Publication, Vol. 388, Submillimetre and Far-Infrared 

Space Instrumentation, ed. E. J. Rolfe & G. Pilbratt, 77 — h 
Crill, B. P., Ade, P. A. R., Artusa, D. R., et al. 2003, ApJS, 148, 527 
Gaertner, S., Benoit, A., Lamarre, J.-M., et al. 1997, A&AS, 126, 151 
Gorski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 
Holmes, W. A., Bock, J. J., Crill, B. P., et al. 2008, Appl. Opt., 47, 

5996 

Huffenberger, K. M., Crill, B. P., Lange, A. E., Gorski, K. M., & 

Lawrence, C. R. 2010, A&A, 510, A58+ 
Jensen, F. & et al. 2010, Proceedings of the Fourth European 

Conference on Antennas and Propagation (EUCAP) 
Jones, W. C, Ade, P. A. R., Bock, J. J., et al. 2006, ApJ, 647, 823 
Jones, W. C, Bhatia, R., Bock, J. J., & Lange, A. E. 2003, 

in Presented at the Society of Photo-Optical Instrumentation 

Engineers (SPIE) Conference, Vol. 4855, Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series, ed. T. G. 

Phillips & J. Zmuidzinas, 227-238 
Lamarre, J., Puget, J., Ade, P. A. R., et al. 2010, A&A, 520, A9+ 
Leahy, J. P., Bersanelli, M., D'Arcangelo, O., et al. 2010, A&A, 520, 

A8+ 

Maffei, B., Noviello, F., Murphy, J. A., et al. 2010, A&A, 520, A12+ 
Mandolesi, N., Bersanelli, M., Butler, R. C, et al. 2010, A&A, 520, 
A3+ 

McDonald, F. B., Webber, W. R., & Reames, D. V. 2010, 
Geophys. Res. Lett., 37, 18101 



20 



The Planck HFI core team: Planck early results: first assessment of the HFI in-flight performance 



Table 6. NEPi is the average of the NEP of processed signal in the band 0.6-2.5 Hz. NEP2 is the white noise component 
of the NEP (see text). The NEATs are derived from the NEP 2 . f is for the bolometers suffering from RTS. 
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